


WEIGHT FUNCTION IN A BIMATERIAL STRIP CONTAINING AN 

INTERFACIAL CRACK AND AN IMPERFECT INTERFACE. 

APPLICATION TO BLOCH-FLOQUET ANALYSIS IN A THIN 

INHOMOGENEOUS STRUCTURE WITH CRACKS. 

A VELLENDER*, G S MISHURIS^, AND A B MOVCHANt 

Abstract. We define a weight function and analyse a problem of anti-plane shear in a bi-material 
strip containing a semi-infinite crack and an imperfect interface. We then present an asymptotic 
algorithm which uses the weight function to evaluate the coefficients in asymptotics of solutions to 

f^^ ' problems of wave propagation in a thin bi-material strip containing a periodic array of cracks situated 

fvQ ' at the interface between two materials. 
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1. Introduction. In this paper we address the problem of determining a weight 
r^ ' function in a domain representing a bi-material strip containing a semi-infinite in- 

O |. terfacial crack. Where the crack is not present the interface is considered imperfect^ 

modelling a thin layer of adhesive between the materials. 

Weight functions are mainly used to evaluate stress intensity factors for asymp- 
totic representations near non-regular boundaries such as crack tips. Classically, sym- 
metric weight functions for interfacial cracks in two-dimensional elasticity were studied 
by Hutchinson et. al [TT] and Bueckner [5]. In these classical works, weight functions 
were defined as the stress intensity factors corresponding to the point force loads ap- 
K*" ' plied to the faces of the crack. More recently, Willis and Movchan [22] defined general 

^^ . weight functions as non-trivial singular solutions of a boundary value problem with 

zero tractions on the faces of the crack and unbounded elastic energy. Recently weight 
functions have been used to perform perturbation analysis of the crack front in |20j 
and to evaluate Lazarus-Leblond constants in J19j . These works contain perfect in- 
terfaces which lead to the well-known square root singularity phenomenon [211 HI] ■ 
In the imperfect interface problem considered in the present paper there is no square 
root singularity in stress components and so the weight function instead takes the role 
of aiding in the evaluation of important asymptotic constants which take the place of 
stress intensity factors. 

^\ ' The imperfect interface is a crucial feature of the problem discussed. Accurate 

^ . asymptotic derivations with various interfaces in composite materials (of imperfect 

type among others) for anti-plane shear without the presence of cracks have been 
analysed in [7] [TUl HI] . Such interfaces have been used to model a thin layer con- 
sisting of small cracks in such a way that the cracks do not appear in the analysis 
in [5l [6j 1^ using the phenomenological approach. Cracks in the static regime with 
imperfect interfaces have been studied in [2j [15] , where it is proved that the imper- 
fect interface leads to a different type of singularity near the crack tip than in the 
ideal interface case. Analysis of the perfect interface with cracks under harmonic load 
can be found in [TJ [T7] and recently for a layered composite with cracks in [9j . The 
manuscript |17| considers wave propagation in a thin bi-material strip and discusses 
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the singular behaviour near the crack tip, while iJLj considers a bi-material plane with 
waves propagating perpendicular to the cracks. 

We consider in this paper Mode III deformation and describe the extent of the 
interface's imperfection by a positive parameter denoted n. The problem we study 
here is a singular perturbation problem; taking very small values for k gives a qualita- 
tively significantly different weight function from that derived for the perfect interface 
case in [17] which corresponds to the formulation with k = 0. Moreover, large values 
of K can lead to interesting effects where the boundary layers surrounding different 
crack tips decay slowly so they can no longer be considered as having no influence on 
the Bloch-Floquet conditions. This effect is discussed in [4]; in the analysis presented 
in the present paper we assume that k is not large enough for these effects to come 
into play and later find a condition for this to be the case. Problems regarding cracks 
in domains including imperfect interfaces have been studied in [5] and }16| . but no 
corresponding weight function has previously been constructed. 

Another critical characteristic of the problem is that the strip considered is very 
thin. In addition to the strip itself being very thin, imperfect interfaces are typically 
replaced with an extremely thin layer of a softer bonding material in finite element 
computations (justified for example in [71 [TOl US]). Moreover, singular behaviour 
exists at the crack tips. These points make FEM modelling for particularly thin 
strips extremely difficult or even impossible and motivate the need for the asymptotic 
approach. In this paper we compare the asymptotic model with finite element simu- 
lations only in cases when the strip is not too thin, but stress that the finite element 
methods are unsuitable for the limiting case whereas the asymptotics remain valid. 
The asymptotic method also obtains crucial constants which describe the solution's 
behaviour at the crack tips which are vital for determining whether fracture may 
occur. These important constants would not be attained by finite element methods. 

The plan of the work is as follows. We first formulate the weight function problem 
and use Fourier transform and Wiener-Hopf techniques |18) to obtain the solution. 
Asymptotic analysis enables us to find analytic expressions for all important constants. 
We then present an application of the weight function to the analysis of Bloch-Floquet 
waves in a structure containing a periodic array of cracks and imperfect interfaces. 
This application involves the derivation of junction conditions. Asymptotic theories 
for structures like rods and plates have received much attention throughout the history 
of elasticity theory. For multi-structures however, conditions in engineering practice 
are often formulated on the basis of intuitive physical assumptions [3] . For example, 
the zero order junction conditions for the problem addressed fit with physical intuition. 
It is important to give these conditions a rigorous mathematical footing; moreover, 
higher order junction conditions do not follow such intuition [12] . 

We conclude by presenting a comparison between the perfect interface case studied 
in [T7] and the imperfect interface case presented here. 

2. Weight Function. 

2.1. Formulation of the Problem. The geometry of the strip in which we 
construct the weight function is shown in Figure [5?T] We define our domain 11^ to be 
the union of 11^ and 11^ , where 

n%^ ^{{x,Y):XeR, (-i)^+ire(o,ff,)}, i = i,2. 

n^ corresponds to the material above the cut with shear modulus /^i, while 11^ 
corresponds to the material below the cut with shear modulus /i2. The materials 
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Fig. 2.1. Geometry for the weight function. 



have respective thicknesses Hi and 7?2- A semi-infinite crack with its tip placed at 
the origin occupies X < 0, while the rest of the interface is assumed to be imperfect 
(see (|2.5p in the text below). 

The functions wi and W2 are defined in domains Hi and 112 respectively as solu- 
tions to the Laplace equation 



(2.1) 



V^Wj{X,Y) =0. 



We impose boundary conditions along the horizontal parts of the boundary of Hb 
and on the crack face itself. We denote the components of stress in the out-of-plane 
direction by 

(2.2) aW(X,y):=M,- 



dr 



J 



1,2. 



We assume a zero stress component in the out-of-plane direction along the top and 
bottom of the strip, as well as along the face of the crack itself: 



(2.3) 
(2.4) 



.(1) 



.(2) 



a'y'l^iX, Hi) = 0, a'^'ziX, -H^) = 0, 



Xe 



4^^(x,o+) = o, 4'l(x,o-) = o, X<0. 



Ahead of the cut we impose the imperfect transmission conditions 

(2.5) «^ily=o+ - ^2\y=o_ = ^41 (X,0+), X > 0, 

where k > is a parameter describing the extent of imperfection of the interface. We 
further assume continuity of tractions across the interface between the materials 



(2.6) 



'yz 



(x,o+) = 4'^(x,o-), X>0. 



We seek solutions in the class of functions that decay exponentially as X — > +oo 
and are bounded as X -^ — cxi: 



(2.7) 



0(e 



-7+X 



), X 



-l-oo: 



= Q+0(eT-^), X^ 



-oo, 



where 7± > and Cj are constants to be sought from the analysis. At the vertex of 
the crack, the solution Wj is assumed to be weakly singular, with 



(2.8) 



wi,W2^0{\n\X\), X^O. 



Formally, conditions (|2.1l) - (|2.7p are similar to those in [T7] if we take n — Q. How- 
ever, with K > the problem is a singular perturbation problem and the behaviour 
described in (|2.8p is entirely different. 
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2.2. An auxiliary problem. We now introduce an auxiliary solution y which 
satisfies the Laplace equation (|2.ip along with the boundary and transmission con- 
ditions (I2.3l) - (|2.6p . but the conditions at infinity and at the vertex of the crack are 
modified as follows: 



(2.9) 

(2.10) 

(2.11) 



= -7+^^, 



(2.12) 



y,^y,io+,o) + o{xin\x\), x^o. 

The functions w and y are related via the formula 

d 



w{X,Y) 



dX 



y{x,Y). 



Bearing this relationship in mind, we often later refer to 3^ as a 'weight function' as 
well as w. It is also shown in [13 that as i? — ;> 0, 



(2.13) 3^± 



i-iya, 



iy) 





TTjJ.j 



/ilKTT 



1 






1-ln 



R 



Rcos9±{7TT0)RsmO} , 



where 3^+ and 3^_ represent yi{R,9) and 3^2 (^,^) respectively and {R,0) describes 
the usual polar co-ordinate system with 9 e [0, tt] for J^i and 9 e [— tt, 0] for 3^2- 

2.2.1. Derivation of Wiener-Hopf equation. We define the Fourier trans- 
forms of yj by 



(2.14) 



3>,(C,y)= j e^^''y,{x,Y)dx. 



The functions y^ are analytic in the strip S = {^ £ 
a double pole only at the point ^ = 0, so 



-7+ < Im(^) < 0}, and have 



yj{s,.Y)^^c,-i^ + o{i), e^o. 



(2.15) 

Note that the functions y, (f , Y) can be analytically extended to the strip 



5 = {?eC:-7+<Im(0<7-}- 
Let us now introduce [3^], the jump in 3^, defined by 

(2.16) [y] = yi\Y^^+-y2\Y^^-- 



We see from (|2.15p that the Fourier transform of the jump [3^](X) generally speaking 
has a double pole at the point ^ = 0. 
We introduce the following notation: 



(2.17)$-(e) = [3^]-Mi«^^ 



Y=0+ 






where we have taken into account (|2.5|) or equivalently the fact that [3^] — Mi^^J- |v-n+ ' 
for X > Q. The function $^(C) is analytic in the half plane Im(^) < and has 
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a double pole at ^ = 0. Thus it can be analytically extended into the half-plane 
C^ == {^ G C : Im(^) < 7_}. We further define the function 



(2.18) $+(^)=/,,|^ 






and so according to (|2.4p . $^(0 is analytic in the half plane C+ = {^ G C : Im(^) > 
-7+}- 



We expect that 



(2.19) <I'^(0 = ^+ ^ r^^ +0(^), e^oo, 






in the respective domain according to (j2.1ip : we later confirm this to be true. 
The Fourier transforms of the functions yj are of the form 

(2.20) y,{^, Y) = AjiO cosh(er) + B,{0 sinh(er). 

Upon the application of boundary and transmission conditions expressions relating 
Aj{£_) and Bj{£_) are found: 

(2.21) B,{0 = {-iyA,{OtanH^Hj), j = l,2; ^iiB,{0 - ^^262(0 = 0. 

Moreover, $^('^) can be expressed in terms of Aj(^), Bj{£_). 

(2.22) $-(o = Ai(e)-A2(o-/ii<si(e), <i>+(o = Mi^i?i(e). 

By applying boundary and transmission conditions, we conclude that the functions 
$+(^) and $~(C) satisfy the functional equation of the Wiener-Hopf type 

(2.23) <i'-(o = -s(e)$+(e), 

in the strip —7+ < Im(^) < 0, where 

(2.24) S(0 = 71— coth(^i/i) + — coth(ei/2) + < 

and —7+ is equal to the size of the imaginary part of the first zero of 5(^) lying below 
the real axis. We would like to stress that the form of the Wiener-Hopf kernel S(^) 
demonstrates that the weight function problem is a singular perturbation problem 
as K -^ 0; the presence of the term involving k fundamentally alters the asymptotic 
behaviour of S(^) as ^ — >■ c». 

2.2.2. Factorization of the Wiener-Hopf kernel. We note that the kernel 
function E!(,^) as defined in (|2.24p can be written in the form 



(2.25) S(0 = . ^^"'^;r "^ ^*(0, 



(A + z^A-zO_ 

e 



where 

^(^1 COth(^iJ2) + /i2 coth(^i7i) + ^1^2kO 



(2.26) S,(e) 



^il^^2Ki\^ + e) 
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and 



(2.27) A^ /Miffi+M2i/2 



Now, S,(^) is analytic in a strip containing the real axis, clearly positive, even and 
smooth for all ^ G R and has been chosen in such a way so that S* (f ) tends towards 
1 as ^ — ?> ±oo and as ^ — > 0. Furthermore, the function S,(^) can be factorized in the 
form 

(2.28) S,(0-S+(OS:(0, 

where 



(2.29) sm-e.^i- J ^^dt 

and /3 > is chosen to be sufficiently small so the contours of integration lie within the 
strip of analyticity of S*(^). The functions SJ^ are analytic in their respective half- 
planes. To conclude this subsection, we have factorised S(^) in the form given in p.25p 
and ()2.28p . where Sj are analytic in the half-planes denoted by their superscripts. 
Note that in the case Hi — H2, other factorisation has been obtained in [2]. 

2.2.3. Asymptotic behaviour of S+. We now seek asymptotic estimates of 
■=■^(0. We first note that for ^ within the strip of analyticity. 



^ + 0(1), S,(e) = l + 0(|en, e^O, ,;=^ + ^ 
? Ml.ni M2W2 



Let us now consider more accurately the behaviour of ^*(^) for ^ e M as ^ — >■ cx). 
Noting that S* (f ) is an even function, it follows from (I2.24p that 



(2.31) S40 = l+ .. .Z -72+0[t7\^], e^ioo. 



The same estimate is true for any ^ lying in the strip of analyticity. We further find 
that 

rrF 

(2.32) EtiO = l + — + Oi\e), C^O, 

m 

-+/fA_i , 1 (Mi +M2)ln(-iO , ^ /^ 1 



7^^ mn2K ^ VKI/ 

the derivation of these expressions is given in Appendix A. Here we have defined the 
asymptotic constant 

00 

(2.34) a = J i^dt. 



The important expression (J2.33I) describing logarithmic asymptotics at infinity is 
needed later for equation (|2.43l) . 
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2.2.4. Solutionof the Wiener-Hopf equation. The factorized equation ()2.23[) 
is of the foriTL 

(2.35) - «(A - zO'^+imtio - Y^e'i'-io^- 

Both sides of (j2.35p represent analytic functions in the strip —7+ < Ini(^) < 7_. 
Moreover we now have asymptotic estimates for 5J^(^) at the zero point in equation 
(|2.32p and for ^ -^ ±00 in (|2.33p . We deduce that since both sides of (|2.35l) exhibit the 
same behaviour at infinity in their respective domains according to (j2.19p . both sides 
must be equal to a constant, which we denote A. We can therefore obtain explicit 
expressions for $^, which are as follows: 

(2.36) *.K) = -^4=^^. «-K) = -^'^+g'--«' . 



We deduce that 

(2.37) yAi,y)- ^^^ I sinHa-iy^'H,) /' '-^'^- 

This allows us to investigate the behaviour of yj as ^ — >■ ±00 and at the zero point. 
It also enables us to find the hitherto unknown real constants Cj and Dj. 

2.2.5. Evaluation of constants Cj, Dj, aj, , 7±. In this subsection we eval- 
uate the constants 7+ (defined in p.9p ). 7_, Cj, Dj (defined in p.lOl) ') and Cg 
(defined in ()2.13p ). We see from our expressions for yj and 4>+ (equations (|2.36p and 
(|2.37p ). along with our asymptotic estimate for S+(^) as ^ ^ that 

i-iy+'^A ( }__l (_a _ 1 
kXii.jHj \e ^\ TT A 



(2.38) y,-(0- \..:. ^ ( 72-7 1 -r-T))+Q(l)^ e^O, 



where a is the constant defined in (|2.34p . It follows from our definition of Cj and Dj 
in (|2T5)) that 

(2.39) C, = 41^, D,^^-:^(^- + \ 

KAfijHj kX/ijHj \7r A 

For normalisation we choose A — kX, giving 

(2.40) C, ^ ^—^ . D, = ^-^ - + - 

The chosen normalisation leaves (|2.40l) i in the same form as in jTTj, but it is clearly 
seen that the expression for Dj (which depends upon k is different). Mishuris (2001) 
[13] demonstrates that near the crack tip (i.e. as i? — > 0), yj{R,d) has behaviour 
described by (12.131) . From this we see that 

(2.41) [y]r^-Ka^^\ R^O. 

The imperfect transmission conditions (j2.5p therefore give that 



(2-42) m ^ 



-4^\ x~^o. 



Y=0+ 
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We earlier made an assumption in (j2.19p regarding the behaviour of $+ at infinity 
and now verify that this was correct. It follows from the expression for $^(^) given 
in (|2.36l) and the asymptotic estimate for S+(^) given in (j2.33p that 



,+ r^N_ ^ , (A'i+M2)A,^, ^^, , ^/^ 1 



(2.43) $+(0 = ^ + ^"^ ,, ln(-»0 + 0{j^], lm(0^+oo, 

which justifies our previous claim. Theorem IB. II (using /^i^ in place of '/' in the 
statement of the theorem) then yields that 

(2.44) ^lim^;.,|^ = -A, 
where A has been defined in (|2.27p and so it follows that 

(2.45) 4^' = A. 

The constant 7+ is the distance of the first zero of 5(^) below the real axis. Manipu- 
lation of (|2.24p indicates that zeros of S(^) satisfy 

(2.46) — cot(7+ffi) + — cot(7+i72) - K7+ = 0, 

Ml fJ'2 

For the first zero below the axis, for large k, 7+ should be small, and so it can be 
shown that 

(2.47) 7+(k) = A(k)(1 + 0(k-1)), k -^ 00, 
indicating that 7+(k) = 0(k^^/^), k^ 00. We also see that 



(2.48) 7+(0) e 

The constant 7_ is given by 



\2Hi' 2H2 



(2.49) ^_=7rmin|— ,— 

In conjunction with (J2.40I) we have now found all constants describing the asymptotic 
behaviour of the weight function 3^. 

3. Application to Analysis of Bloch-Floquet Waves. In this section, we 
present an application of the weight function derived in the previous section by ad- 
dressing the problem of out-of plane shear Bloch-Floquet waves within a thin bi- 
material strip containing a periodic array of longitudinal cracks and imperfect inter- 
faces. The problem addressed is an imperfect interface analogue to that studied in 

3.1. Geometry. The geometry of an elementary cell of the thin periodic struc- 
ture considered is shown in Figure [3Tl The elementary cell is of length a and contains 
two materials of thicknesses eHi and eH2 , where £ is a small dimensionless param- 
eter. These materials occupy respective domains 11^^', j = 1,2, and the elementary 
cell is further split into smaller domains flj^ , m = 1,2,3,4, as shown in Figure IXTl 
Along the join of the two materials and centered on the origin sits a crack of length 
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n 



(1) 



n, 



(2) 



Fig. 3.1. Geometry of the elementary cell. 



I. Outside the crack, the interface is assumed to be imperfect, which models a thin 
layer of adhesive joining the materials together. The extent of this imperfection is 
represented by the parameter k. 

The functions u'-^^{x, y) are defined in He , j = 1, 2 as solutions of the Helmholtz 
equations 



(3.1) 



V^u'^^\x,y) + '^u'^'\x,y)=Q, {x,y) ^n[^\ j = l,2. 



Here, Cj = \/ lijj Pj are the shear speeds in their respective domains j = 1,2. The 
functions u^^> are regarded as out-of-plane displacements, \Xj denotes the shear modu- 

lus and pj the mass density of the material occupying n^ . The quantity w represents 
the radian frequency of the time-harmonic vibrations with amplitude u. 

3.2. Boundary conditions. We impose boundary conditions along the horizon- 
tal parts of the boundary of Hg and on the crack face itself. We use similar notation 
to that in the previous section to denote the components of stress (see (|2.2p ). 

We assume a zero stress component in the out-of-plane direction along the top 
and bottom of the strip, as well as along the face of the crack itself: 



(3.2) 
(3.3) 



4i)(x,eiJi) = 0, .,, 



(2) 



aW(x,0^ 



0, 



(x, -eiJa) =0, xe (-a/2, a/2), 
ag)(x,0-)=0, a;G(-//2,//2). 



Outside the crack, along the boundary between He and He , there is an imper- 
fect interface described by the condition 

(3.4) u(i)(x,0+)-u(2)(x,0~) =eKCTW(a;,0+), xG (-a/2, -^2) U (//2, a/2). 
We also assume continuity of stress across the interface 

(3.5) aW(a;,0+) = af;(x,0-), a; G (-a/2, -//2) U (Z/2, a/2). 

We seek the solutions u^^' which represent the Bloch-Floquet waves, so that at the 
ends of our elementary cell x — ±a/2 we have for j = 1,2 the Bloch-Floquet conditions 



(3.6) 

(3.7) 



7/W(-a/2,y) 



_ p~iKa Aj) 



rU)l 



u^^>ia/2,y), 



(-a/2,y) = e-'^V«(a/2,2/), 



ye {-eH2,eHi), 
ye {-eH2,eHi). 
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For a fixed value of tlie Bloch parameter K, we seek the eigenvalues w and the 
corresponding eigenfunctions u'^'^^ with finite norm in W2(J^' ), j — 1,2. 

In p.4p . the case in which k = corresponds to an ideal/perfect interface between 
the different materials; such a problem was considered in [T7]. Where possible we will 
follow the same line as in this paper. To summarise the approach, we approximate u 
in a certain form, derive a lower-dimensional model together with boundary layers in 
the vicinity of the vertices of the crack and then use our weight function to assist in 
the derivation of junction conditions for a skeleton model. 

3.3. Asymptotic Ansatz. The eigenfunctions u{x, y) are approximated in the 
form 

k=0 lm=l 

(3.8) + (W^^^ {Xa,Y) + wI,^\Xb,Y)^} + RN{x,y,e), 

with scaled co-ordinates Xa, Xb and Y introduced in the vicinity of the left and right 
vertices of the crack defined as 

(3.9) Xa^^—^, Xb^^-^^, Y=y. 

e £ e 

(rn) . . . .... (7) 

Here, Wj, represent solutions of lower-dimensional problems within limit sets ilj, , 
j = 1, 2, 3, 4. Xm = Xm{x, y, e) are cut-off functions defined so that Xmix, y;e) = l in 
f2e and decay rapidly to zero outside fJe . They vanish near the so-called junction 
points A and B (the vertices of the crack) . The terms W]^ and Wg represent the 
boundary layers near A and B, and Vm is the 'fast' change of eigenfunctions in the 
transverse direction in the domain fl^ . R^ is the remainder term in the asymptotic 
approximation. We would like to indicate to the reader that the uppercase scaled 
co-ordinate Xb defined in p.9p corresponds to X from the derivation of the weight 
function in section [5] 

We note that this form of Ansatz relies upon the vital assumption that the bound- 
ary layers surrounding the crack vertices A and B are independent. That is, we assume 
that the exponential decay of both boundary layers is sufficiently rapid so that it is 
negligible in the vicinity of the other crack tip. 

In this paper we will consider the form of approximation given in p.Sp with A^ = 1 
and will comment on the effect of taking higher order approximations. 

3.4. One-dimensional model problems. Outside the vicinity of A and B, 
the boundary layers W)^ and Wg decay (we later verify this to be the case) and so 
seek u in the form 

(3.10) u{x,y,e) ^ ^ e'^ (t;W(x) + ^^^^(a;, F)) , 

fe=0 

where V,n have zero average over the cross-section of fi^ ^ for all ttt, = 1,2, 3, 4. That 

is, 

Hi ° 

(3.11) J \i^Hx,Y)dY = 0, J Vi^Hx,Y)dY = 0. 

-H2 
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Since the low-dimensional model problem studied in 17 was the same above and 
below the crack (in f2e ' and 51^ '), we refer the reader to that paper. The problem 
is however differently formulated in f2e and f2e due to the imperfect transmission 
conditions in these domains. We focus our attention on the layered structure Qe , 
analogous arguments will apply to i7e . We use the notation i)J ■ to denote the 
function v^ in He . The key observation is then to note that the transmission 



condition across the imperfect interface as given in (I3.4p imply that 
(3.12) wj^"^ _ w(^) = 0, fc = 0, 1. 

and so it follows that for k = 0, 1, that the solution to this low dimensional model is 
not impacted by the presence of the imperfect interface. 

To conclude this section, we have found that our case with the imperfect interface 
has the same equations for the low dimensional model up to terms in e as the case 
with the perfect interface studied in [IT]- The equations for vl and V^ are of 
course similar to the case examined here where m — 1. We would like to stress that 
the imperfect interface impacts on the low dimensional model equations for terms 
in e'^, k > 2. The equations gained in this section need to be complemented with 
the boundary conditions and junction conditions at the points xa and xb- In order 
to derive these junction conditions which depend on the imperfect parameter k, we 
construct boundary layers in the vicinity of the vertices of the crack. 

4. Junction conditions. We introduce four smooth cut-off functions Xm € 
C°°(M) in the spirit of [17]. These are functions defined so that Xm{x,y;s) = 1 in 
fie' and decay rapidly to zero outside fl^- These allow us to extend the function 
([XTUl) outside rji™', m = 1, 2, 3, 4, giving 

(4.1) u(x,y;£) ~ ^ ^' E Xm{x,y,e) [v^^Hx) + e^Vi^Hx,Y)) , 

k=0 m=l 

however this gives an error near the junction points xa and xb- We therefore introduce 
boundary layers Wa{Xa,Y) and Wb{Xb,Y), and so seek u{x,y,e) in the form 

X™ (u(? (x) + e^V^^ {x, Y)) + W^''> {Xa, Y) + w'i^ {Xb,Y) \ . 

Substituting this expression into the original equation and comparing terms of the 
same degree of e we obtain 

(4.2) V^^y{w^W(X„,y)+J-W(X„,r)}=0, a = A,B, A; -0,1, 
with the functions Fa ,fc = 0,l, a = A, i? given by 

3 4 

^T = E '"^n{xA)Xm{x,y;£), J"jj°' == ^ v^^\xB)Xm{x,y\e), 

771—1 771 — 2 

3 

(4.3) T^l^ = Y, [{v^^^)\xA)XA+vil\xA)}xm{x,y;e), 

7n— 1 
4 

(4.4) 4^) = Y. {{v^^^)'{xB)XB+v(^\xB)]xm{x,y;e). 
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Fig. 4.1. Contour of integration for p 



We now focus our attention near xb] analogous arguments apply to xa- We will con- 
sider in the following analysis four functions gi, i — 1,2, 3, 4, which are solutions of the 
Laplace equation. These solutions also satisfy the boundary conditions corresponding 
to zero stress on the top and bottom edges of the strip (|2.3p as well as along the cut 
itself (j2.4p . They also satisfy the transmission condition (|2.5I) across the imperfect 
interface, along with continuity of stress (j2.6p . These solutions are given by 



(4.5) 



51 = 1, 52 = Xi 



93 = y, 



9i = 



dy 



where y is the weight function derived in section [51 



(fc) 



Since they are boundary layers, we expect that W^ decay exponentially as X 



i(fe) 



D)'"' asX ^ -cx). We first express C)'"' , D)'"' , fc = 0, 1 



(fc) 



+0O and behave as C- 'X 

in terms of «,„ and their derivatives. We have from Green's formula that 



iW i.- 



(4.6) 



= 



2 

E 
=1 



^ij 






4'V(<^+-^. 



-{k).dg. 



on 



an^^'(L) 



The further analysis is quite similar to that in |17| . although we would like to stress 
that the weight function y in the present paper is different, as are the transmission 
conditions. We therefore need to prepare this analysis from the beginning where it is 
different for g^, g^. 



4.1. The cases k = 0,1, 



1,2,3. We see from boundary conditions that 



integrals over the horizontal parts of the boundary l\ ,lr^ ,1-4 



;(j) jU) lU) 



J 



1 , 2 give zero 
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contribution to the integral. Moreover, the contribution from Sg also disappears as 

li^'' and lij'' 



5^-0 (see Figure ICT for gi, g2 and 53, leaving contributions solely from /j ^-nd l^ 



in these cases. 

From the definitions of J^^ , we obtain the following limits at ±(X) for k = 0, 1: 

(4.7) tI°^ = vf\xB), Xb -^ +00, 

(4.8) ^ij"^ = 4°) (xB)H(r) + 4"^ (a;s)H(-y), Xs ^ -00, 

(4.9) J-jj'^ = (wr)'(a;B)^B + v[^\xb), Xb^+oo, 

3 

(4.10) 4')=^{(«f)'(a;B)XB+«f^(xB)}H((-l)^r), Xb^-oo, 



j=2 



(k) 



where HiY) is the Heaviside step function. Since Wg — ?> as Xb — > +00, equation 
(|4.6p reduces to 



'" , 



2 



("^) - E'. / (».4 (^?' + »<•) - (^1=" + < ') ^) <<« 



*9 



Applying this procedure with each of 31,52, <?3 and Tb\^b yields six equations, 



which are presented in subsection 14.31 

4.2. The cases fc = 0,1, i = 4. To obtain a further two equations, we apply 
the same procedure to the solution 54 = ^^ . Again, the contribution from the 
horizontal parts of the contour of integration is zero, leaving nonzero contributions 
from the vertical parts of the contour, l^ and l\^ . Unlike with gi, 52 and 53 however, 
the contribution from Sg^-'' is non-zero. We investigate the behaviour of 54 near the 
crack tip. 

We have that 54 = -^ = "^^ cos0 — -^-^ sinfl, where {R, 9) is the usual polar 
co-ordinate system, with R = \/ X"^ + Y"^ and so from our asymptotic estimate for 
yj near the crack tip from we deduce that near the crack tip, 

(4.13) 5?^ ^ ^^ 14^^ +4^hni?+ (-l)(^+i)4^^in20(^ + (-1)^-0)1 , 

and so for small R, 



(4.14) 



OR TT/iji? 



(i) 
Noting that the outward normal to Sg is in the direction of —R, we have that as 



S^O 






(4.15) 




Ss 


= 


Ss 


Ss 
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Since Wb satisfies the same model problem, as y, it too will possess asymptotic 
behaviour at the crack tip of the same form as 54 in (|4.13p , but with different constants 
which we denote aLJ and ^L-i for fc = 0, 1. The contribution to the integral from 
the circular part of the contour is therefore given by 



jy) 



(y) 



(wl,'\o+,e) + ^|j''(o+, 0)) de + ^ (wi'\o+,0) + :fP{o+, 0)) de 



Ay) 



^0 



TT J TT^l 1 + ^ e^) 




Ay) 



— 7r 




M2 






(y) {w) 



B ' 



With this information at hand, we are now able to apply (|4.6p with 54 and J^g \j'\ 
yielding a further two relationships. 

4.3. Deriving the junction conditions. We define the column matrices 

T 



(4.16) 



i^(fc) 



(j{k) (j{k) 



Df^ D. 



(k) 



0,1. 



The eight equations obtained in the previous two subsections can then be rewritten 
as two matrix equations, the first of which is found to be 



(4.17) M^(°) 





(^ii?i + lJ.2H2)vf' {x b) - iiiHivf\xB) - tJ^2H2vf\xB) 

lilH^Civf\xB) + ^i2H2C2vf\xB) 

^jy)..iw) 



where M is the 4x4 matrix 



(4.18) 



HiHi 


/^2^2 














fiiHi 


A'2^2 


lnHiDi 


mH2D2 


-fiiHiCi 


— /^2-'^2C2 



filHiCi /i2^2C'2 











where Cj and Dj are the asymptotic constants from the weight function defined in 
dSini)- The determinant of M is given by det(M) = -fil^j,lH^H^{Ci - 02)'^ < 0. 
Therefore for c|°^ = cf ^ = D^"^ 



D2 = (that is, for W to vanish far away from 



the crack tip as we would expect for such a boundary layer), we have that the matrix 
in the right hand side of (|4.17l) must be equal to zero. From this follow the junction 
conditions 



(4.19) 
(4.20) 



(0)/ N {0)f \ {0)f N 

^2 i^B) = ^3 [XBJ = ^4 [XBJ, 



AW) 
'■(0) 



0. 



The latter condition (|4:20l) yields that W^°^ 



(4.21) 



B 



0. The second matrix equation is 
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{HiHi + ^J.2H2)v[^ (xb) - /iii?iW2 (a^s) - ^2^2^ (^^s) 

^ilHlClvi'\xB) + mH2C2vi^\xB) - fiiHiDi{v!,°'^y{xB) - fi2H2D2{v^°^Y{xB) 

Ka^^'^a^^^ - fi,H,Ci{vi''^y{xB) - M2i?2C2(t^f )'(xb) 

where M is the matrix given in (|4l8l) . For c|^^ = C;^^^ = D^^^ = Z)^^^ = 0, the 
right hand matrix is again set to zero. Noting that a^ ' = A (see (|2.45p ) and that 
fiiHiCi + fJ,2H2C2 = 0, setting the fourth row of the RHS matrix to zero then yields 
that 

(4.22) a;^' = ;^A{(.(0))'}. 
where 

(4.23) A{(«(o))'}(xs) = ivi'^YixB) - (t'f )'(^b). 
The other conditions imply 

M2^2 fa V 

(4.25) 4"(«) - .;■'(«) + jj^, {I + i) A(("'"')'}(«). 
along with the relationship 

(4.26) {iJiiH^+ ^l2H2){vf'^y{xB) - ^llH^{v^°^y{xB) - ^l2H2{vf^y{xB) = 0. 

We stress that a and A are functions of k and so expressions (|4.24l) and (I4.25P de- 
scribe how the junction conditions depend upon the extent of imperfection of the 
interface. In particular, {aj-n + 1/A) is a constant that plays a crucial physical role 
since it defines the proportionality between the displacement jump in the first order 
approximation and the angle of opening in the zero order approximation. Equation 
(|4.26p complements conditions (|4.19p and (I4.20p to give full information for the zero 
order approximation. We later present numerical results for the normalized constant 
ai^[a/^) + l/\)l{Hi+H2). 

The conditions regarding the first order approximation (|4.22p . (|4.24p and (|4.25p 

can be complemented by a further equation in (wm )'(2:s), which follows from the 
next level of approximation, i.e. taking iV = 2 in p.Sp : 



(4.24) v^^\xb) = v['\xb) .. /"". ^ ( 3 + 4 ) A{(«(°))'}(XB) 



2^2 

^ilH,{v^^^yixB)+^i2H2{v^3^y{xB)~{fllH^+^l2H2){v^^^yixB)^Y. / ^^B^^n^'^ 



and by our earlier comment that W^ — 0, the right side of this expression is zero. 
At this point we would like to comment that taking higher order approximations and 
evaluating higher order junction conditions is possible but much more advanced. For 
example, integrals analogous to that on the right hand side of the above expression 
would depend upon Wj^ and boundary layers from higher order approximations, and 
so would not in general be zero. However, since we focus on thin strips, e is small 
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Fig. 5.1. Contour plots of the ratio aj/ap for four different values of «;«, a dimensionless 
parameter describing the extent of imperfection of the interface between the two materials. The 
axes of each plot are fit, and Ht, dimensionless parameters respectively describing the mechanical 
and geometric properties of the problem. The ratio aj /ap gets closer to 1 as k, decreases in value 
towards 0. 



and so terms in e^ would give significantly less contribution than the lower order 
approximations. We later comment on the accuracy of the zero order approximation 
on page [19] by comparing computations against FEM results in a case where e is not 
too small. We would like to underline that the accuracy will increase for smaller e, 
but for very small e it is no longer possible to obtain finite element computations. 

5. Numerical simulations and discussions. To enable us to compare results 
with the perfect interface case discussed in [T^ effectively, we seek normalized con- 
stants. We first seek a normalized representation of a. We introduce the notation 



H = H,+ H,. H. = |L-#, ^. . ^^l^i^, .. = <tl±M^ A. = XH, 



H1+H2 



Ml + M2 



H 



where i/* , /i* and k* are non-dimensional parameters which respectively describe the 
geometrical, mechanical and imperfect properties of the problem. A is the constant 
dependent on /i^, Hj and k defined in (J2.27I) . A* can be expressed in terms of the 
other dimensionless parameters as 



(5.1) 



xi 



8{1 + ^i,,H^) 
«.(1-m2)(1-H2)- 



We also introduce the function 
t ( 2 



^■ifi-^ \tj 



t + 



XI +t'^ V ny + ^J^.*) 



coth 



t{l+H,) 



hi:^ [ i fl:^ j 



coth 



i(l - ff*) 
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which satisfies the relationship S^H,(i) — S, (-4), and so we can write 
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(5.2) a = \'^,i- 





^ -H / dt ~ Ha^,, 



t2 



H 



t2 



where we have defined the non-dimensional quantity a^. We find through asymptotic 
analysis that 



(5.3) 



ln5„(i) 1 iJ,V* 



H? 



M*ff* + 1 



i2 



12 



1 + fj,i,H^ 



0{t^), t^O. 



Mishuris, Movchan and Bercial [T7| showed that in the analogous problem to that 
discussed in this paper with a perfect interface instead of an imperfect interface, 



(5.4) 
where 



ap = — In 

TT 



D,^ap{Hi+H2)Cj, 



1 + H^ 1-H^ 

l + HA ■'- fl-H^' 






if* — tanh(t7:f*) coth(t) 
(sinh(t) + /i, sinh(ti/»))i 



dt. 



We have demonstrated (see the form of the constants Cj, Dj in (j2.40p ) that for the 
imperfect interface problem, 



(5.5) 
where 



D^^ai{Hi+H2)C, 



J' 



(5.6) 



1 f \nE,,{t) ^^ 1 



and since small k* correspond to an interface which is 'almost perfect', we would 
expect a I — > ap as k* — >■ 0. Figure [STTj shows a plot of the ratio aj/ap on axes of /i* 
against i?, for four different values of k*. From this it is easily seen that as k* — >■ 0, 
ai/ap gets close to 1 as expected. The behaviour of the weight functions near the 
crack tip are however absolutely different since the problem is singularly perturbed, 
that is: 



(5.7) 



[y] 



frjR, K — ?► 0, 



(5.8) 



Ml 



dy_ 

dY 



F=0+ 



K 



where -q is defined in (|2.30p . In the plots showing the ratio for small values of k*, the 
highest deviations from 1 occur near the corner of the plot. These correspond to the 
cases where there is a high contrast between the shear moduli and thicknesses of the 
two materials. We see that in the cases where the materials have similar shear moduli 
and thicknesses (nearer the center of the plot), the ratio ai /ap quickly approaches 1 
as K» — >■ 0. 
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Perfect case ( ap) 





Fig. 5.2. Surface plots of af for k, = 100, 1 and 0.01; also of ap, all plotted on axes of fi^ 
and Hf . 



Figure [5T2] shows surface plots of aj on axes of /i* and if* for k* = 100, 1, and 
0.01. This constant describes the impact that the imperfect interface has upon the 
junction conditions as described in equations (|4.24l) and (|4.25p . Also shown in the 
figure is a plot of ap. The similarity between the plot of aj for k* = 0.01 and the 
plot of ap is evident here. For the cases with larger k* values, we see that ap is 
differently dependent upon the mechanical and geometric parameters of the problem. 

Figure [531 shows finite-element plots (COMSOL) of standing wave eigensolutions. 
For these simulations we use the following geometrical parameters for the elementary 
cell: 



I = 0.8[m], 



2.4[m], Hi = 0.1[m], H2 = 0.05 [m], 



t(2)^ 



and the following material constants which correspond to iron (in He ) and aluminium 

(inni^'). 

Ai2 = 82 • 10^[N/m2], ^1 = 26 • lO^iN/m^], p2 = 7860[kg/m3], pi = 2700[kg/m3]. 

Presented in this figure are three plots corresponding to Al-Fe strips with different 
materials bonding them together, with the vertical dotted lines indicating the location 
of the crack tips. The imperfect interface is modelled in the COMSOL simulations by 
a thin layer occupied by an adhesive material; this approach was justified in |15l 116] . 
among others. Provided that h^csin/ H2 is sufficiently small and /^rcsin is small in 
comparison to fii and /i2, this gives k — /ircsin/z^rcsin- 
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Fig. 5.3. Finite element computation (COMSOL) contour plot of the eigensolution correspond- 
ing to the standing Bloch- Floquet waves for three different values of k. Top: Bonding material with 
shear modulus lOOO/^i-csin • Middle: Bonding material is epoxy resin. Bottom: Bonding material 
with shear modulus ^rcsin/10. Countours join points of integer values, and the dotted vertical lines 
indicate the location of the crack tips. 



The second of the three plots in Figure 15.31 uses epoxy resin as the bonding 
material with parameters 



A^r 



= 2.5-10^[N/m2 



= 1850[kg/m3], /i,e,in = 0.01[m]. 



For comparison, the first plot shows a simulation with a gluing layer of shear modulus 
1000 greater than that of epoxy resin. The third plot uses a material with shear 
modulus 10 times less than epoxy resin. Equivalently, these three cases in the top, 
middle and bottom parts of the figure correspond to k» = 2.88 • 10^'^, k, = 2.88, and 
K* = 28.8, respectively. The plots show that the standing wave is more localised and 
intense in the locality of the crack when the bonding material is stiffer. Conversely, 
when the bonding material is less stiff, the standing wave extends further beyond the 
locality of the crack and is less intense. Closely packed contours indicate areas where 
stress is high; we see that the highest stress is to be found in the vicinity of the crack 
tip in all three cases. Moreover, as we would expect, the highest stress intensity is 
found in the case with the stiffest bonding material. 

We do not present dispersion diagrams here computed by the asymptotic analysis 
and COMSOL as they are similar to those given in [Ij. As in that paper, the biggest 
discrepancy between results obtained from asymptotic analysis and numerical simula- 
tions appear for the case of the standing waves. In all other situations the accuracy is 
very good, with a typical discrepancy between finite element and asymptotic results of 
around 0.3% in the case where the strip has the same dimensions as used throughout 
this section, which corresponds to e = 0.0625. We remind the reader that we use 
static boundary layers in the analysis. The standing waves lie in the area of rather 
high frequencies, which may provide one possible explanation for this phenomenon. 
This discrepancy needs to be eliminated and this will form part of future work. 

It is readily seen in the bottom plot of Figure 15.31 (which corresponds to a highly 
imperfect interface) that the boundary layer support extends almost to the edge of 
the elementary cell. This extension far away from the crack tips suggests that the 
boundary layers decay slowly from the crack tips and so may not be assumed indepen- 
dent. In this case, therefore, our analysis may become invalid due to the assumption 
in our asymptotic procedure that the exponentially decaying boundary layer does not 
influence the Bloch-Floquet conditions. This assumption is satisfied if 7+ is far from 
zero, so if k is not too large. More accurately, we assume 7+ ^ ^^ (see (|2.47p for 
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large k). If the imperfect interface is too weak and this condition is violated then 
the junction conditions evaluated here will no longer be accurate and other analysis 
should be sought. 
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Appendix A. Derivation of asymptotics of S+(^). We present here the 
derivation of asymptotics for S;^(^). The results of this derivation are used in expres- 
sions (|2.32|) and ()2.33|) . We introduce the auxiliary function 

oo— i/3 



(A.i) e+(e)- ] -jz^dt, 

— oo—i(3 

so that S+(^) = exp((l/27ri)e+(C)) (see ((2:29)) ). We first note that e+(0) = since 
the integrand is odd and estimate (|2. 301) 9 demonstrates integrability of S* at the zero 
point, allowing us to take /3 = 0. Thus 

oo — iP oo — ip 



e+(0 = 

— oo — z/3 



lnS,(f) lnS,(i) 



since the integral is bounded. Also, we have that 



dt^^ I ^^df^O, C^O, 



(A.2) / '^dt^J'-^dt = 2j'^dt^2a, 

— oo — ifi — oo 

since the integrand is even and again by considering (J2.30p 9. which indicates that we 
have integrability at the zero point. Here we have found that 

(A.3) e+(c) = 2< + o(|ep), e^o. 

From this we obtain the following estimate for S+(^) as ^ — >■ 0: 

(A.4) s+(o = i+^+o(i^n, e^o. 



TTl 



We now seek estimates of 6j(C) for ^ — >■ oo within the domain. To avoid problems 
caused by integrating along the real line, we consider ^ -^ cx) in such a way that 
Im(^) — ^ +00. Integrating (jA.ip by parts, splitting the integral in two and manipu- 
lating the resulting expression gives 





We introduce an arbitrary i? > and split this integral at R to give 

R 

We then see that 

(A.7) ^"(t^)='^ + ^(|f)' ^^^' ^<'<^^ 
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and from (j2.3ip we have 

Kit) 



This aUows us to estimate 



A^i +M2 



K^lfl2 



^ + o(^], t^oo. 



QtiO 






'■"l^)^' + Kl'- «-°°' 



After integrating by parts and performing a change of variables, we find that 



(A.9)|ihr 



1 , fi+1 



dt 



1 


/ 


1 




In 




c 


V 


e 



«arg( -Tj) ) +0( — ) , C^oo, 



and so from (|X|, we deduce that 
2(mi +Ai2) 



(A.io) e:(0 



hrHC)+0 - , lm(0^+cx.. 

MiM2< VkL 



Recalhng the relationship between our auxiliary function 9;,^ and S+ as we discussed 
after (lA.ll) . we see that 

(A.ll) ^+(0 = 1 + — T + 777 , lm(0^+oo. 



Appendix B. Theorem. 

Theorem B . 1 . Let f(x) be the function 



(B.l) 



fix) = ^ I <i>+(t)e-"*di. 

Zn 



If^^{t) is analytic in C^ and 

(B.2) $+(f)=ait-i+0(f-(i+^)), t^oo, 

where S > is small, in the closed half-plane C = C^ U M, i/ien /(a;) = for all 
X < and 



(B.3) 



lim f(x) = —iai. 



Proof. The fact that f{x) = for all a; < is a direct consequence of the fact 
that <&^(i) is a '+' function. Assume now that 2: > 0. From the assumptions on 
the behaviour of the function $^(t), it follows that ^~^{t) = ait~^ + R{t), where 
tR{t) -^ 0, as t -> 00, t e C (including t -^ ±cx), t e M). 

We write 



(B.4) fix) 



{00 a ^ 

J J 

a —a J 
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The first integral is 

CX) 

(B.5) /i(x, a) = A$+(-t)e'-* + $+(i)e-"*]di = fii{x, a) + /i2(x, a), 

a 

where 

(B.6) fii{x,a)-- 



' t t 



, sin(i) , 
dt = -2iai I — ^-^di, 



and 



/i2(x,a) = / [i?(-t)e"* + i?(t)e-"*] di 



t 



X \ X I X \x _ 



dt. 



Taking a — a;~^/^, we have that fii{x,x^^^^) -^ —inai and fi2{x,x~^/^) — >■ as 
a: — >■ 0+. Let us denote the second integral in (|B.4[) by f2ix, a). Then using analyticity 
of ^'^{t) in C+ and defining 

(B.7) Ta^{te C\t ^ ae'^ < 6* < tt}, 

we deduce 

(B.8) f2ix,a)^- U+{t)e-"'dt. 

We write this in the form 

(B.9) f2{x,a) ^ f2i{x,a)+ f22ix,a), 

where 



(B.IO) f2iix,a) 



f ^e-"*di, and/22 (x, a) = - I R{t)e-~'^'dt. 



Again taking a = x ^'^, we obtain 
{B.n)f2i{x,x-^'^) = - J ^e-"*dt: 



-ai 



—dt = —inai, x — > 



-1/2 



r -1 



/2 



Now, 
(B.12) 



f22(x,a) 



f R{t)e-''''dt = ~ I -R 



1 „ ft 



e~'*dt, 



and so f22ix,x ^^'^) ^- as a; — >■ 0+. By collecting together these observations and 
reconsidering equation (|B.4|) we conclude that 



1 



(B.13) 

which completes the proof. D 



f{x) — > -— (— iTTfli — iTTai) = —iai, x — >■ 0^, 



